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In this work, we developed a numerical analysis-associated experiment method to determine the effective mul- 
tiplication factor kett, which is difficult to obtain directly from conventional neutron source multiplication (NSM) 
method. The method is based on the relationship between ker, subcritical multiplication factor ks and external 
neutron source efficiency * in the subcritical system. On basis of the theoretical analysis, the dependence of 
ks and &* on subcriticality and source position was investigated at the Chinese Fast Burst Reactor-II (CFBR-II). 
A series of k, were measured by NSM experiments at four subcritical states (Ker = 0.996, 0.994, 0.991 and 
0.986) with the ?5?Cf neutron source located at different positions (from the system center to outside) at each 
subcritical states. The * was obtained by Monte-Carlo simulation for each condition. With the measured ks 
and calculated ®*, keg of the subcritical system was evaluated with a relative difference of < 1% between values 
obtained by the improved method and by positive period method. Especially, the relative difference of < 0.18% 


with the source located at the system center. 
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I. INTRODUCTION 


Neutron source multiplication (NSM) method, a simple 
measurement technique that uses an external neutron source 
and an ordinary neutron detector, is suitable for real-time 
measurement. With the recent development of accelerator- 
driven system (ADS) [1], an increasing attention has been 
paid to NSM experiments associated with various subcritical 
systems [2, 3]. Implementing this kind of experiments can ac- 
quire subcritical multiplication factor ks and external neutron 
source efficiency &*, which can be used to design the ADS 
systems, and its neutron multiplication [4, 5]. The subcritical 
multiplication factor k, is related to the neutron multiplica- 
tion M [6] by ks = 1 — 1/M, but in practice one needs to 
know the effective multiplication factor keff- 

In this paper, in order to overcome the difficulties of getting 
kett directly from the conventional neutron source multiplica- 
tion method, we propose a numerically associated experiment 
method. In Sec.II, the NSM method and the relationship be- 
tween effective multiplication factor keg, subcritical multipli- 
cation factor k; and neutron source efficiency &* are reviewed 
briefly. The experiment apparatus is described in Sec. III. In 
Sec. IV, for verifying the proposed NSM method, the numer- 
ical analysis and a series of NSM experiments were carried 
out at the Chinese Fast Burst Reactor-II (CFBR-ID at various 
subcritical states with ?°?Cf neutron source located at differ- 
ent positions in each subcritical states (from center of the sys- 
tem to outside). 


* Supported by Natural Science Foundation of China (No.11175164) 
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Il. THEORETICAL BASIS 


For the critical systems composed of the fissile materials 
and other materials, the neutron balance equation in funda- 
mental mode is expressed as: 


1 


Lyo(r, E) = t a 


Pyo(r, E), (1) 


where, yo(r, Æ) is the neutron flux distribution at position 
r(x, y, z) at energy E, L is the destruction operator, and P is 
the production operator. By taking the inner product < ... > 
of Eq. (1), which indicates the integration over the space and 
energy variables, kere in Eq. (1) can be expressed as: 


_ < Pyo(r, E) > 


kefe = ; 2 
ee Lyo(r, E) > 2) 


For the subcritical system with an external source at steady 
state, the neutron balance equation is: 


Ly, (r, E) = Pys (r, E) + s(r, E). (3) 


Similarly, the subcritical multiplication factor k, can be ex- 
pressed as: 


< Pys(r, E) > 

he Se N (4) 
< Lys(r, E) > 

where, y.(r, E) is the neutron flux in the subcritical system 

with the external neutron source s(r, E). 

The external source efficiency, &*, is the ratio of average 
importance between the source neutron and fission neutron. 
The &* implies that a neutron emitted from source, with a 
certain position and energy distribution, will produce certain 
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“important” neutrons that can lead to higher or lower neutron 
multiplication for a given proliferation system. 

Deduced from properties of the operator and Eqs. (1)- (4), 
we have Eq. (5): 


ks 


k- Fk- © 


1 : @* (1 : k 
keff 7 ( x) ve: 
Then, the multiplication factors k; and keff can be theoreti- 

cally correlated with the neutron source efficiency &*. 
Obtaining kere of the system through Eq. (5) needs to know 

the k, and 4*. As mentioned above, the k; can be mea- 
sured by neutron multiplication experiment, but the quantity 
of neutron source efficiency &* cannot be obtained from ex- 
periments directly. In the following sections, we acquired &* 
from numerical analysis and evaluated ker successfully. 


II. EXPERIMENT APPARATUS REVIEW 


Experiments were carried out at the Chinese Fast Burst 
Reactor-II (CFBR-ID [7], an ellipsoidal assembly, which con- 
sists of, from inside to outside, the enriched-uranium spher- 
ical shells, inner brass reflector, depleted uranium reflector 
and outer brass reflector. The reactor can be divided into 
the upper, middle and lower parts (Fig. 1). The middle part 
is a 5.2cm thick stainless steel plate, which has three hori- 
zontal holes for three control rods to pass through, i.e. an 
auto- adjustment rod, a composition rod and a pulse rod. 
The three rods are made of casted enriched-uranium. A 
10cm x 4.2cm x 14.5cm (w x h x d) irradiation chamber 
and two horizontal experimental holes are located in the plate. 
The neutron source can be put at any position in the horizon- 
tal holes by a neutron source transmission device. Reactivity 
of the assembly can be controlled by adjusting insert depth 
of the three control rods. The assembly works at subcritical, 
delayed critical or super prompt critical states. 

A *°2Cf spontaneous fission neutron source (in activ- 
ity of 1.41 x 10~°s~+) sealed in an aluminum box of 
P7.8 mm x 11mm can be sent into a horizontal hole of the 
stainless steel plate by the transmission device when the sys- 
tem is in deep subcritical state. Then, the control rods are 
inserted to different depths to operate the reactor at different 
subcritical power states. 

Two BF3 proportional counters, of the SZJ-1 type, pro- 
duced by Beijing Nuclear Instrument Factory, are used to 
record the neutrons (which is proportional to the neutron den- 
sity in the system). They consist of the paraffin barrel (Han- 
son long counter), preamplifier, HV power supply, main am- 
plifier and multi-scalar. The detectors are 180 cm away from 
the system center. 


IV. NUMERICAL ANALYSIS AND EXPERIMENTS 
A. Numerical calculations 


Before a series of experiments are carried out, Monte-Carlo 
numerical method [8] was used to build the subcritical sys- 
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Fig. 1. Schematic of CFBR-IL. 


tem model that is identical to CFBR-II in Fig. 1. Four sub- 
critical states of ket =0.996, 0.994, 0.991 and 0.986 were 
selected, which are conditions for our multiplication experi- 
ments. The source positions were changed along the r-axis 
at each subcritical state (i.e. 0mm, 20mm, 40mm, 60mm, 
80mm, 100mm, 120mm and 150mm). The spontaneous 
252Cf neutron source, described by Watt fission spectrum [9], 
is positioned the same way as that in the planned experiments. 
The ks of all the experimental conditions were calculated by 
MCNP5 with statistic difference of < 0.7%. The neutron 
source efficiency &* is closely connected with neutron impor- 
tance, which is proportional to the detector response to a unit 
source. In other words, neutron importance is a measure of 
the “importance” of a neutron in contributing to the detector 
response, or the expected counts per neutron in certain posi- 
tion, direction and energy. The detector counting rates to the 
system with the 7°2Cf neutron source in the core could be ob- 
tained, and the counting rate to the eigen neutron source dis- 
tributed in the system could be recorded, too, by using MCNP 
code. The ratio of the two counting rates should be the ex- 
ternal neutron source-related average efficiency &* [10, 11]. 
In terms of eigen-distribution neutron source acquisition, we 
used the KCODE and SSW cards in MCNPS critical calcu- 
lation to achieve the goal. Thus, we tallied * in MCNP5S 
calculation following certain regulation according to MCNP 
rules. 

By calculating the ks, &* and kett inside and outside of the 
neutron source, we were able to investigate the effect of sys- 
tem subcriticality and source position. 

For the subcritical states (Keg = 0.994), the effect of the 
source positions along the r-axis are demonstrated by the cal- 
culation results in Fig. 2. The source efficiency ®* and sub- 
critical multiplication factor k, at different source positions 
indicate that the smallest discrepancy between k, and keff is 
achieved when the source is at the center of the core. 

Table 1 shows the subcriticality effect on k, and & at 
kett =0.996, 0.994, 0.991 and 0.986. One sees that for the 
four subcritical states, when the source is at the core center, 
ks is very close to ker, with a relative difference of < 0.8%, 
and the &* was greater than 0.91. The results indicate that the 
closer the system reaches to the critical state, the more like- 
liness the relative values of y, and yo become equal to each 
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other, and so do the values of k, and kere. Inversely, placing 
the source outside the core, the discrepancy between k, and 
kett becomes much larger, with the maximum relative differ- 
ence being 22.1340% at kep = 0.986 and the source being 
150 mm away from the core center. 


Kep= 0.994 


kett ks 
© 
n 


P 
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Fig. 2. (Color online) Effect of source position on subcritical mul- 
tiplication factor ks and source efficiency ®* at a fixed subcritical 
state of ke = 0.994. 


TABLE 1. Summary of subcriticality effect on ks and ©* for a ?°?Cf 
neutron source in the subcritical model 


case kett Source position (mm) (ks /kett — 1) (%) p* 
1 0.9960 0 —0.6436 0.9599 
150 — 12.2234 0.0251 
2 0.9940 0 —0.5851 0.9272 
150 — 15.2628 0.0887 
3 0.9910 0 —0.7513 0.9471 
150 — 18.6502 0.0939 
4 0.9860 0 —0.7852 0.9176 
150 —22.1340 0.1036 


In order to check the source position effects on the neutron 
leakage rate, we did another numerical analysis to investi- 
gate the detector response to non-collided neutrons emitted 
from the source, and to total neutrons. The results in Fig. 3 
show that the relative difference is below 0.25%. Placing the 
neutron source closer to the system edge, the counting rate 
for non-collided neutrons from the source becomes greater, 
while the counts of total neutrons decrease. In other words, 
the detector counting rate decreases for the eigen-distribution 
neutron (fission neutron inside the system). This explains 
why the discrepancy between k, and keff becomes large while 
moving the ?5?Cf neutron source away from the system cen- 
ter. 
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Fig. 3. (Color online) The detector counts for total and non-collided 
neutrons (kett = 0.994). 
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Fig. 4. (Color online) ks measured at different subcritical states. 


B. Experiments Results and Discussion 


On the basis of the numerical investigation in Sec.IV.A, the 
NSM experiments were carried out at CFBR-II to measure ks 
and evaluate *. The subcritical states were experimentally 
achieved by adjusting insert depth of the control rods. All the 
experiment apparatus and conditions were completely identi- 
cal to the numerical ones. The results are shown in Fig. 4. The 
ks decreased with increasing distance of the source from the 
core center and with deepening subcriticality. This is identi- 
cal to the calculation results in Fig. 2. 


In order to test the feasibility of this improved NSM 
method, we considered four subcritical states (Ker = 0.996, 
0.994, 0.991 and 0.986) and eight neutron source positions 
according to each subcritical states. For the sake of the calcu- 
lated source efficiency * and the measured ks, we obtained 
the effective multiplication factor keg according to Eq. (5) at 
32 studied conditions. For comparison, we used the positive 
periodic method to measure the reactivity as the reference 
method in the work. The specific process is that we measured 
reactivity in super-criticality, and extrapolated to subcritical- 
ity with calibrated reactivity curve of control rods. 
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Fig. 5. (Color online) Experiments and calculations results (superscript “C”: Calculated value; superscript “E”: Experiments value; superscript 


“ke: Evaluated value). 


As shown in Table 2, a comparison between the reference 
method and the results in this work at the four subcritical 
conditions with the source located at the core center revealed 
fairly good agreement with a maximum relative difference of 
0.18%. For a certain subcriticality (Ker = 0.991), the discrep- 
ancy between kett measured by positive period method and 
the evaluated ones in this work is the smallest of all. When 
kett > 0.9910, the evaluated kett are a bit larger than experi- 
mental measured ones. However, at a very deep subcritical- 
ity, e.g. Keg = 0.9860, the evaluated kert is a little less than 
the measured one. This is because the source neutrons could 
not propagate efficiently in deep subcriticality. The results 
demonstrated the good precision of this improved method and 
the validation of the Monte-Carlo calculation model. 


In general, kere of the system for all studied cases were well 
approximated as shown in Fig. 5. 


The comparison between the Keg results in refer- 
ences (positive period method, marked as keff—exp) and 
our method (keff—eva) is done by Ae = [(keff-eva — 
Kett—exp)/Kett—exp] xX 100%. As shown in Fig. 6, the space 
effect on keff can be seen in a large range of the neutron posi- 
tion. Generally, when the neutron source is inside the subcrit- 
ical system, the relative difference is small, comparing to that 
when the source is outside the system; and the relative dif- 


TABLE 2. Comparison between measured and calculated kett, ks and 
®* in the fast critical assembly experiments while neutron source put 
at the center of core. 


Case kete? keff k p* ( kerf” Ikeff > 1) 1% 
1 0.9960 0.9977 0.9936 0.9599 0.1707 
2 0.9940 0.9948 0.9899 0.9272 0.0805 
3 0.9910 0.9915 0.9847 0.9471 0.0505 
4 0.9860 0.9844 0.9757 0.9176 —0.1623 


* Positive periodic method. 

> Calculated in this work. 

© Neutron Source Multiplication experiments and Eq. (1). 
4 Numerical analyzed method described in Sec.IV.A. 


ference increase rapidly with distance of the neutron source 
to the system center, being irregularly large at the distance of 
> 85cm, because this is where the brass and depleted ura- 
nium reflectors locate, causing uncertainty factor related to 
complicated neutron behavior. 


Therefore, to reduce the relative difference in experiments 
implemented, the neutron source should be put inside the sub- 
critical system. 
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Fig. 6. (Color online) Relative difference at four subcritical states 
between calculated and measured ones. 
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V. CONCLUSION 


A methodology to evaluate the effective multiplication 
factor keg from measured k, and numerical analyzed ®* 
was investigated via NSM experiments on CFBR-II. The 
evaluated results were approximately equal to the reference 
method (positive period method). 

On basis of the theoretical and numerical preparations, ks 
was measured through a series of NSM experiments with a 
252Cf spontaneous fission neutron source located at different 
positions in four subcritical states of CFBR-II. The effective 
neutron multiplication factor ker obtained in all the subcriti- 
cal states when neutron source located inside of the subcrit- 
ical system, within a relative difference of < 0.2% between 
the values obtained by the improved NSM method and pos- 
itive period method. So, to reduce the relative difference in 
experiments implemented, the neutron source should be put 
inside the system. The numerical analysis-associated NSM 
we proposed is feasible for the simple subcritical apparatus, 
such as CFBR-II. 
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